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Abstract 

The neutrino oscillation experiment KamLAND has provided us with the first evidence for v e 
disappearance, coming from nuclear reactors. We have combined their data with all solar neutrino 
data, assuming two flavor neutrino mixing, and obtained allowed parameter regions which are 
compatible with the so-called large mixing angle MSW solution to the solar neutrino problem. 
The allowed regions in the plane of mixing angle and mass squared difference are now split into 
two islands at 99% C.L. We have speculated how these two islands can be distinguished in the near 
future. We have shown that a 50% reduction of the error on SNO neutral-current measurement 
can be important in establishing in each of these islands the true values of these parameters lie. We 
also have simulated KamLAND positron energy spectrum after 1 year of data taking, assuming the 
current best fitted values of the oscillation parameters, combined it the with current solar neutrino 
data and showed how these two split islands can be modified. 
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I. INTRODUCTION 

Many solar and atmospheric neutrino experiments have collected data in the last decades, 
giving evidence that neutrinos produced in the Sun and in the Earth's atmosphere suffer 
flavor conversion. While the atmospheric neutrino results [l] may be understood by — > v T 
conversion driven by a neutrino mass squared difference within the experimental reach of 
the accelerator based neutrino oscillation experiment K2K [2], the mass squared difference 
needed to explain the solar neutrino data was, until quite recently, before the Kamioka 

n 

Liquid scintillator AntiNeutrino Detector (KamLAND) |3| has started its operation, too 

small to be inspected by a terrestrial neutrino oscillation experiment. 

A number of different fits, assuming standard neutrino oscillations induced by mass and 
1 fl 
mixing as well as other exotic flavor conversion mechanisms |5| , have beenperformed using 

the combined solar neutrino data from Homestake fl, GALLEX/GNO Q Q, SAGE Q, 

Super-Kamiokande-I [l(3] and SNO jllj]. These analyses selected some allowed areas in the 

free parameter region of each investigated mechanism, but did not allow one to establish 

beyond reasonable doubt which is the mechanism and what are the values of the parameters 

that are responsible for solar flavor conversion. After the first result of the KamLAND 

(or KL hereafter) experiment 3( this picture has changed drastically. 

In the first part of this paper, we present the allowed region for the oscillation parameters 

in two generations for the entire set of solar neutrino data, for KamLAND data alone and 

for KamLAND result combined with all solar neutrino data, showing that this last result 

finally establishes the so called large mixing angle (LMA) Mikheyev-Smirnov-Wolfenstein 

(MSW) ^| solution as the final answer to the long standing solar neutrino problem jl^ . 

definitely discarding all the other mass induced or more exotic solutions. (For the first dis- 

n 

cussions on the complete "MSW triangle" which includes the LMA region, see Ref. [lj.) 
In the second part, we speculate on the possibility of further constraining the oscillation 
parameters in the near future. For instance, we point out the importance of SNO neutral- 
current (NC) data in further constraining the LMA MSW solution. In particular, we discuss 
the consequence of a significant reduction (50 %) of the SNO neutral-current data uncer- 
tainty. Finally, we simulate the expected inverse /3-decay e + energy spectrum after 1 year 
of KamLAND data taking, based on the best fitted values of the oscillation parameters. 
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We combine this with the current solar neutrino data in order to show how the allowed 
parameter region can be modified. 



II. DETERMINATION OF OSCILLATION PARAMETERS 

KamLAND has observed about 40% suppression of v e flux with respect to the theoretically 
expected one |3|, which is compatible with neutrino oscillations in vacuum in two generations. 
In this case the relevant oscillation parameters, which must be determined by the fit to 
experimental data, are a mass squared difference (Am 2 ) and a mixing angle iff). We first 
obtained the allowed region in the (tan 2 8, Am 2 ) plane compatible with all solar neutrino 
experimental data, then with KamLAND data alone, and finally we combine these two sets 
of data. 



A. Solar Neutrino Experiments 



We have determined the parameter region allowed by the solar neutrino rates measured by 
Homestake fl, GALLEX/GNO Jjj, SAGE B 

and SNO (elastic scattering, charged-current 
and neutral-current reactions) llj (6 data points) as well as by the Super-Kamiokande-I 
zenith spectrum data (44 data points), assuming neutrino oscillations in two generations. 

We have computed the v e — > v e survival probability, properly taking into account the 
neutrino production distributions in the Sun according to the Standard Solar Model |l5|. 
the zenith-angle exposure of each experiment, as well as the Earth matter effect as in Ref. |5|, 
except that here we solved the neutrino evolution equation entirely numerically. We then 
have estimated the allowed parameter region by minimizing the x% function which is defined 
as 



Xq 



(J. 



pth no 

R j -Rj 



obs 



(1) 



M=l,...,50 

where Rf 1 and R° hs denote the theoretically expected and observed event rates, respectively, 
which run through all 50 data points mentioned above, and a & is the 50 x 50 correlated 
error matrix, defined in a similar way as in Ref. In this work we have treated the 8 B 
neutrino flux as a free parameter. 
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In Figure Q we show the region, in the (tan 2 9, Am 2 ) plane, allowed by the Super- 
Kamiokande-I zenith spectrum data as well as by the rates of all other solar neutrino ex- 
periments at 90%, 95%, 99% and 99.73% C.L. In our fit we obtained a ^(min) = 37.7 for 
47 d.o.f (83 % C.L.), corresponding to the global best fit values Am 2 = 7.5 x 10~ 5 eV 2 and 
tan 2 6 = 0.42. 

B. KamLAND 

KamLAND [3] is a reactor neutrino oscillation experiment searching for v e oscillation 
from over 16 power reactors in Japan and South Korea, mostly located at distances that 
vary from 80 to 344 km from the Kamioka mine, allowing KamLAND to probe the LMA 
MSW neutrino oscillation solution to the solar neutrino problem. 

The KamLAND detector consists of about 1 kton of liquid scintillator surrounded by 
photomultiplier tubes that register the arrival of v e through the inverse /3-decay reaction 
v e +P —> e + + n, by measuring e + and the 2.2 MeV 7-ray from neutron capture of a proton 
in delayed coincidence. The e + annihilate in the detector, producing the total visible energy 
E which is related to the incoming v e energy, E u , as E = E v — (m n — m p ) + m e , where m n , 
m p and m e are respectively, the neutron, proton and electron mass. 

After 145.1 days of data taking, which corresponds to 162 ton yr exposure, KamLAND has 
measured 54 inverse /?-decay events, where 87 were expected without neutrino conversion. 
These events are distributed in 13 bins of 0.425 MeV above the analysis threshold of 2.6 
MeV (applied to contain the background under about 1 event). 

We have theoretically computed the expected number of events in the z-th bin, N^ heo , as 

Af co = / dE v a{E u )Y,4>k{Ev)P^v e [ dE R(E,E'), (2) 

k Ji 

where R(E, E') is the energy resolution function, E the observed and E' the true e + energy, 
with the energy resolution 7.5%/y^(MeV). Here cr(E u ) is the neutrino interaction cross- 
section and 0fc is the neutrino flux from the k-th power reactor, we have included all reactors 
with baseline smaller than 350 km in the sum. P Ve ^ Ve = Pp e ^p e (if CPT is conserved, which 
we will assume here) is the familiar neutrino survival probability in vacuum (the matter effect 
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is negligible here), which is equal to one in case of no oscillation, and explicitly depends on 
Am 2 and tan 2 9. 

We were able to compute the region, in the (tan 2 6, Am 2 ) plane, allowed by the Kam- 
LAND spectrum data, by minimizing with respect to these free parameters, the Xkl function 
defined as Xkl = Xg + Xp with 

XI = J2 { 1 / l ] , (3) 

and 

jyobs 

Xp = E 2(iVf e ° - Nf) + 2 iVf s In (4) 

3 3 



where = y N° bs + (0.0642 N° bs ) 2 is the statistical plus systematic uncertainty in the 
number of events in the i-th bin and the sum in is done over the bins having 4 or more 
(less than 4) events. We have also computed the allowed regions using purely Gaussian 
or Poissonian \ 2 functions and found that the hybrid \ 2 definition above could reproduce 
better KamLAND's allowed regions jjj. Therefore, we have prefered to use it in our paper 
(see also Ref. [lfijl). 

Using this Xkl we have computed the allowed region at 90%, 95%, 99% and 99.73% C.L. 
shown in Fig. 12 which are quite consistent with the ones obtained by the KamLAND group 
in Fig. 6 of Ref. Q]. In our fit we obtained a XKL( m i n ) = 5-4 for 11 d.o.f (91 % C.L.), 
corresponding to the best fit values Am 2 = 7.0 x 10~ 5 eV 2 and tan 2 6 = 0.79. 



C. Combined Results 

Combining the results of all solar experiments with KamLAND data we have obtained 
the allowed region showed in Fig. El The minimum value of xiot = X% + Xkl f° r ^ ne 
combined fit is Xtot( mm ) = 43.6 for 60 d.o.f (94.5 % C.L.), corresponding to the best fit 
values Am 2 = 7.1 x 10~ 5 eV 2 and tan 2 # = 0.42. We observe that there are two separated 
regions which are allowed at 99 % C.L. : a lower one in Am 2 (region 1) where the global 
best fit point is located, and an upper one (region 2) where the local best fit values are 
Am 2 = 1.5 x 10~ 4 eV 2 and tan 2 9 = 0.41, corresponding to Xkjc( m i n ) = 49.2. We observe 
that depending on the definition of Xkl (gaussian, poisson or hybrid) used, a third tiny 
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region above Am 2 = 2 x 10~ 4 eV 2 appears at 99.73% C.L. However, apart from this small 
change, the combined allowed region is not essentially affected by the Xkl used. 

In Fig. |3] we show the theoretically predicted energy spectra at KamLAND for no oscil- 
lation, the best fit values of the oscillation parameters for KamLAND data alone and for 
KamLAND combined with solar data in regions 1 and 2. We note that the fourth energy 
bin, which is for the moment below the analysis cut, can be quite important in determining 
the values of the oscillation parameters in the future. 

III. FUTURE PERSPECTIVES 

In this section we consider the effect of possible experimental improvements which can 
help in determining the oscillation parameters with more accuracy in the future. We first 
consider a reduction of the error in the SNO neutral-current measurement then an increase 
of event statistics in KamLAND. 

A. Effect of reducing SNO neutral-current error 

In order to constrain the solar neutrino oscillation parameters even more, in particular, 
to decide in which of the 99% C.L. islands Am 2 really lie, we have investigated the effect 
of increasing the SNO neutral-current data precision to twice its current value. We have 
re-calculated the region, in the (tan 2 9, Am 2 ) plane, allowed by all current solar neutrino 
data, artificially decreasing the SNO NC measurement error but keeping the current central 
value, as well as the other solar neutrino data, unchanged. The result can be seen in Fig. |SJ 
The best fit point and the value of ^©(niin) remain practically unchanged with respect to the 
result obtained in Sec. Ill A[ but the allowed region shrinks significantly. This is because the 
8 B neutrino flux normalization, which can be directly inferred from SNO NC measurements, 
gets more constrained. Combining this with KamLAND data we obtain the allowed region 
shown in Fig. El We observe that this allowed region is substantially smaller compared to 
the one shown in Fig. |3J Moreover, region 2 only remains at 99% C.L. 
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B. Effect of increasing KamLAND statistics 

We simulate the expected KamLAND spectrum after one year of data taking for three 
distinct assumptions. We have generated KamLAND future data compatible with the best 
fitted values of Am 2 and tan 2 9 obtained for : (a) KamLAND data alone, (b) KamLAND 
and current solar neutrino data in region 1 and (c) KamLAND and current solar neutrino 
data in region 2. We have also included an extra bin, corresponding to the fourth bin in 
Fig. HJ We have re-calculated the region allowed by the combined fit with the current solar 
neutrino data in each case. 

The results of our calculations can be seen in Figs. [70 If the future KamLAND result 
is close to the current one (see Fig. EJ), values of tan 2 # larger than the ones allowed now 
will be possible and region 2 will be excluded at 99% CL. For this case we have obtained 
X 2 ot (min) = 42.1. On the other hand, if the future KamLAND data are more compatible 
with the current best fit point of solar neutrino data (see Fig. [HJ, the global allowed region 
will diminish substantially with respect to Fig. El and region 2 will only remain at 99% 
C.L. For this case we have obtained x 2 ot (min) = 39.1. Finally, if after one year KamLAND 
data is more compatible with region 2 (see Fig. then one should observe an increase 
towards larger values of Am 2 in the combined allowed region with respect to the one shown 
in Fig. El In this case region 1 and 2 will have similar statistical significance, corresponding 
to Xt 2 ot( min ) ~ 44 -0- 



IV. DISCUSSIONS AND CONCLUSION 

We have performed a combined analysis of the complete set of solar neutrino data with 
the recent KamLAND result in a two neutrino flavor oscillation scheme. We have obtained, 



in agreement with other groups |l7j . two distinct islands, denominated regions 1 and 2, 
in the (tan 2 9, Am 2 ) plane, which are the most probable regions where the true values of 
these parameters lie. Region 1, where the global best fit point was found, is around Am 2 = 
7.1 x 1CT 5 eV 2 , while region 2 is around Am 2 = 1.5 x 1CT 4 eV 2 . 

We have considered two possible future improvements in the determination of the neutrino 
oscillation parameters. First, we have investigated the effect of a 50 % decrease in the error 
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of the SNO NC measurement. We have shown that this would substantially reduce the 
allowed parameter region when combined with KamLAND data. In particular, region 2 
would not be allowed at 95% CL. anymore. 

Second, we have studied what can happen in the near future, when KamLAND collects 
1 year of data. We have simulated the expected KamLAND spectrum including an extra 
lower bin, corresponding to the fourth bin in Fig. EH Three different cases were studied in 
combination with the present solar neutrino data. In the first case, we have assumed that 
the future KamLAND spectrum will be compatible with oscillation parameter values at the 
best fit point for the present KamLAND data alone. This is the most restrictive case for 
region 2. In the second case, we have considered that future data will be more compatible 
with the present best fit point for the solar neutrino experiments. In this case, the combined 
allowed region will be much smaller than the present one and region 2 will be only allowed 
at 99% CL. Finally, in the third case, we have assumed that the future KamLAND data will 
be compatible with the local best fit point in region 2. In this case, the combined allowed 
region will suffer an increase towards larger values of Am 2 and region 1 and 2 will both have 
similar statistical significance. 
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FIG. 1: Region in (tan 2 9, Am 2 ) plane allowed by the Super-Kamiokande-I zenith spectrum com- 
bined with rates from Homestake, GALLEX/GNO, SAGE and SNO. The best fit point is marked 
by a cross. 
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FIG. 2: Regions in (tan 2 6, Am 2 ) plane allowed by KamLAND data alone. The best fit point is 
marked by a cross. 
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FIG. 3: Region allowed by all the solar neutrino experiments combined with KamLAND (KL) 
data. The region below (above) Am 2 = 10 -4 eV 2 is referred to as region 1 (2). The best fit points 
in each region are also marked by cross (global best) and plus (local best). 
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FIG. 4: Expected positron energy spectra at KamLAND (KL) for no oscillation, the best fit values 
of the oscillation parameters for KamLAND data alone and KamLAND data combined with the 
solar neutrino data in regions 1 and 2 of Fig. |31 The KamLAND data 3( is a ls° shown as solid 
circles with error bars. The energy threshold at 2.6 MeV is marked by a vertical line. 
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FIG. 5: Same as Fig. ^ but decreasing the SNO neutral-current data error to half of its current 
value. 
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FIG. 6: Same as Fig. |2l but decreasing the SNO neutral-current data error to half of its current 



value. 
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FIG. 7: Same as Fig. El but for a simulated KamLAND spectrum after one year of data taking 
compatible with the KamLAND alone best fit Am 2 = 7 x 1CT 5 eV 2 and tan 2 9 = 0.79. 
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FIG. 8 
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Same as Fig. [7|but for the KL + Solar neutrino global best fit Am 2 = 7.1 x 10 5 eV 2 and 
= 0.42 in region 1. 
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t: Same as Fig. [7| but for the KL + Solar neutrino local best fit Am 2 = 1.5 x 10 4 eV 2 and 
= 0.41 in region 2. 
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